A MODEL DESIGNED TO PREDICT THE 
MOTION OF OBJECTS TRANSLATED 
BY CLASSICAL BLAST WAVES 


|. Gerald Bowen, Ray W. Albright, 
E. Royce Fletcher, and Clayton S. White 


Issuance Date: June 29, 1961 


CIVIL EFFECTS TEST OPERATIONS 
U.S. ATOMIC ENERGY COMMISSION 


= 


NOTICE 


_ This report is published in the interest of providing information which may prove of 
value to the reader in his study of effects data derived principally from nuclear weapons 
tests, | s oe 
This document is based on information available at the time of preparation which 
may have subsequently been expanded and re-evaluated, Also, in preparing this report 
for publication, some classified material may have been removed. Users are cautioned 
to avoid interpretations and conclusions based on unknown or incomplete data. 


-” PRINTED IN USA | | ‘ 
Price $1.25. Available from the Office of 
Technical Services, Department of Commerce, 
Washington 25, D. C. 


USAEC Office of Technical Information Extension 
Oak Ridge, Tennessee 


A MODEL DESIGNED TO PREDICT THE 
MOTION OF OBJECTS TRANSLATED 
BY CLASSICAL BLAST WAVES 


By 


I. Gerald Bowen, Ray W. Albright, 
E. Royce Fletcher, and Clayton S. White 


Approved by: R. L. CORSBIE 
Director 
Civil Effects Test Operations 


Lovelace Foundation for Medical Education and Research 
Albuquerque, New Mexico 
January 1961 


co) 


ny 


ABSTRACT 


A theoretical model was developed for the purpose of predicting the motion of objects 
translated by winds associated with “classical” blast waves produced by explosions. Among 
the factors omitted from the model for the sake of simplicity were gravity and the friction that 
may occur between the displaced object and the surface upon which it initially rested. Numeri- 
cal solutions were obtained (up to the time when maximum missile velocity occurs) in terms 
of dimensionless quantities to facilitate application to specific blast situations. The results 
were computed within arbitrarily chosen limits for blast waves with shock strengths from 
0.068 to 1.7 atm (1 to 25 psi at sea level) for displaced objects with aerodynamic characteris- 
tics ranging from those of a human being to those of 10-mg stones and for weapon yields at 
least as small as 1 kt or as large as 20 Mt. 
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Chapter | 
INTRODUCTION 


1.1 OBJECTIVES 


During the 1955 and 1957 Test Operations at the Nevada Test Site, the masses and veloci- 
ties of over 20,000 objects (window-glass fragments, stones, spheres, sticks, etc.) which were 
translated by nuclear-produced blast waves were experimentally determined,’ along with a 
time-displacement history of an anthropometric dummy simulating man.’ The availability of 
such a mass of data stimulated an analytical study calculated to arrive at a mathematical for- 
mulation capable of predicting the translation of objects by blast, particularly since this data 
offered an empirical fabric against which to test the success of the analytical effort. 

The purpose of this report is to describe, step by step, the theoretical studies that have 
resulted in a mathematical model capable of predicting the motion of objects utilizing selected 
basic blast parameters. This model, however, is applicable only to those situations in which 
*classical wave forms exist. 


1.2 SCOPE AND LIMITATIONS 


The applicability of the model itself has no well-defined limits; however, the numerical 
solutions that were obtained and are reported herein have been arbitrarily limited in scope. In 
general, the aim was to compute velocity, displacement, and acceleration as a function of time 
for objects ranging in size from a pea to man; these computations were to cover biast waves 
with shock overpressures from 1 to 25 psi (14.7-psi ambient pressure) and weapon yields from 
1 kt to 20 Mt. 

Another class of limitations is invoked not by the scope of the computations, but by the 
model itself. Formulation of a workable model was facilitated by the use of certain simplify- 
ing assumptions. These assumptions, which are discussed below, have not, in general, caused 
serious discrepancies between predicted velocities and those measured in the field operations, 
particularly in those situations where the blast wave was classical.t 

As a practical approach, it was assumed first that the effect of surface friction was neg- 
ligible. It has been observed that fairly large objects tend to be lofted when subjected to blast 
waves; the more intense the blast, the heavier the object that it is capable of lifting against 
gravity. Nonspherical objects could develop either positive or negative lift depending on their 
orientation to the wind. Thus, the validity of the no-friction assumption is dependent upon the 
strength of the blast wave, the object under consideration, its random orientation, and the na- 
ture of the surface over which translation occurs. It will be shown later that certain uses can 
be made of the data even for situations in which surface friction is a significant factor. 


*The term “classical blast wave” is used in this report to mean the typical wave not ap- 
preciably modified by terrain effects and possessing a well-defined shock front. 

+A limited discussion of the agreement between predicted and measured velocities is, 
made later in this report. A more complete treatment will be found in Ref. 3. 
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A second approximation made concerned the assumption that there was no gain or loss of 
energy as a result of the object’s moving with or against gravity. The kinetic energy that is 
lost during lofting would be regained as the object fell to its original elevation, thus mitigating 
somewhat the error in the predicted motion. 

Third, only the propelling force of the wind was considered. Another force that might have 
been included was that due to differences in overpressure between one side of the object and 
the other during passage of the shock front (diffractive loading). Since the bodies being con- 
sidered were relatively small (up to the size of man), the classical blast wave would engulf the 
object very quickly and impart only a small momentum as a result of the overpressure itself. 

Fourth, it was assumed that there was no change in the properties of an object which gov- 
erned acceleration (area presented to the wind, drag coefficient, and mass) during the accel- 
erative phase of displacement. For irregular, rigid objects that are nearly spherical, such as 
stones, this is a reasonable assumption. For objects that are obviously nonspherical or de- 
formable, prediction of a range of velocities taking into account both maximum and minimum 
drag areas is often used. Another useful procedure is to employ the average drag area de- 
rived from the concept of random orientation.® 

Fifth, no allowance was made for the fact that a displaced object may be moved to a lower 
overpressure region and thus be acted upon by correspondingly weaker blast winds. The re- 
sults of the computations themselves seem to justify the neglection of this phenomenon. It 
will be shown that displaced objects receive a large percentage of their velocities in a rela- 
tively short distance over which the decay of the blast wave is small. 
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Chapter 2 
ANALYTICAL PROCEDURE 


2.1 NOMENCLATURE 


The terminology used in this study is defined in this section. A lower-case letter is used 
to represent a quantity with dimensions. In general, the same letter is capitalized if the quan- 
tity is made dimensionless by an appropriate factor or factors. Thus, the dimensionless term 
is represented by its principle variable. The factors, or parameters, used to make quantities 
dimensionless invariably are constants for any given blast situation. 


a = acceleration coefficient = sCy/m 
A = apoty/co 
Cy = speed of sound in undisturbed air 
Cq = drag coefficient of moving object 
d = distance of travel of moving object 
dm = distance of travel of moving object when maximum velocity is reached 
D = d/(tic,) a 
I = overpressure impulse = {‘? P dt 
m = mass of moving object 
p = overpressure or pressure in excess of py 
Po = pressure of undisturbed air 
Pp, = maximum or shock overpressure 
P=p/py 
Ps =P, /Po; shock overpressure in atmospheres 
q = dynamic pressure = (1/2)pu’ 
qd, = dynamic pressure at the shock front 
Q =q/Ppo 
Q, = 4,/Po 
p = air density 
Ss = area presented to wind by moving object 
= time after arrival of blast wave 
tS = duration of positive pressure phase of blast wave 
t{ = duration of winds in the direction of propagation of the blast wave 


T =t/tp 
u = velocity of the air 
U= u/c 


v = velocity of the moving object 
Vm = Maximum velocity of moving object 


V= v/Co 
v = acceleration of moving object 
v= vti/c, 


W = weapon yield in kilotons 
x = velocity of propagation of the pressure disturbance 
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x = X/Cp 
Z=t/tt 


NOTE: Any variable that is underlined indicates the average value taken over a particular 
time interval. 


2.2 EQUATIONS OF MOTION 


2.2.1 Fundamental Concepts 


Newton’s second law of motion can be stated as 


_. av 
Force =m at 


and the drag force on a moving object is 
1 2 
Force = 5 plu —v)* sCq 


since the net wind moving past the object is (u ~ v). Combining the above equations and solving 
for dv, we obtain 


ay =p —_ at (2.1) 


It was convenient to isolate and label the physical parameters that involve the moving ob- 
ject in Eq. 2.1. 


sCd_» (2.2) 


Now, @ can be called “acceleration coefficient” since it completely describes the object in so 
far as the computation of velocity vs. time is concerned. Thus, two objects possessing the 
same value of a, regardless of dissimilarity of shape, size, and mass, would experience the 
same increase in velocity if exposed to the same or similar blast waves. 


2.2.2 Time Correction 


Since the moving object, or missile, travels along with the blast wave, the time during 
which it is exposed to blast winds is longer the higher its velocity is in relation to the velocity 
of propagation of the pressure disturbance and associated winds. Consider a small segment of 
the blast wave of length dx where the air-particle density and velocity are approximately con- 
stant. If this segment moves with a velocity x, then 


dx = x at 


where dt is the time required for the segment dx to pass a fixed point. Similarly, the velocity 
of propagation of a blast-wave segment past a missile that itself is moving at velocity v is 


(x —v), and 
dx = (x — v) dt’ 


where dt’ is the time required for the blast segment to pass the missile. By eliminating dx 
between the above equations, we obtain 


dt’ = dt — (2.3) 
xX —V 
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Combining Eqs. 2.1 and 2.2 and substituting the corrected time dt’ of Eq. 2.3 for dt, we obtain 


dt (2.4) 


dv =F pa(u —y)? 


x 
(x —v) 


It was more convenient to work with dynamic pressure, q = (1/2) pu’, than with air den- 
sity. For this reason Eq. 2.4 was modified to 


_ u—v\? x at (2.5) 
av = qa ($Y) &—Vv) 


2.2.9 Dimensional Analysis 


The blast-wave variables in Eq. 2.5 are determined as a function of time by four parame- 
ters: (1) shock overpressure, p,; (2) ambient pressure, pp; (3) duration of the positive winds, 
tt; and (4) speed of sound in the undisturbed air, cy. Computations were made for particular 
values of shock overpressure in atmospheres, P, = p,/Po- The last three parameters, how- 
ever, were used to make the variables in Eq. 2.5 dimensionless. The obvious advantage of this 
procedure is that computed values of missile velocity, distance of travel, and acceleration can 
be modified after the computations have been made to fit any blast situation defined by po, t7, 
and Co. 

The variables of Eq. 2.5 were made dimensionless through the application of the following 
algebratic operations: (1) both sides of the equation were divided by cy, (2) the numerators and 
the denominators of the two fractions were divided by cy, (3) @ was multiplied by p,) and q was 
divided by py, and (4) t was divided by tt and @ was multiplied by tj. 

After these operations have been performed, Eq. 2.5 becomes 


a5) = (B) CEES) [ee] [eres tres «(G) a 


and, after appropriate substitutions (see Sec. 2.1, Nomenclature) 
_ u-v\? x 
av = aa( U xy 92 (2.7) 


Two additional quantities are used in dimensionless form, distance of travel and accelera- 
tion. Since both are functions of velocity and time, their dimensionless forms are determined 
by dimensionless velocity and time. Thus 


= — Vv t = de oa 
D=VZ = tt Catt (2.2 
and 
- V_vtt_ vt 
=— ae = 2.9 
7 2 Cot 0 ( ) 


2.2.4 Approximation Solution 


The explicit expressions of Q, U, and X as a function of time for a particular blast wave 
are very cumbersome. Added to this difficulty is the fact that the variable V cannot be sepa- 
rated from the time-dependent variables (see Eq. 2.7). Hope for a complete mathematical so- 
lution was soon abandoned. 

A stepwise solution was then attempted which would permit the blast parameters to be 
held constant for small increments of time but would allow the missile velocity to vary as a 
nonlinear function of time. So that a simple mathematical integration can be accomplished, the 
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time-correction term, X/(X — V), was not included. Thus 


Vg tAV ZytAZ’ 
dv 
a 2 


Vo 


where V,) and Z, are the velocity and time, respectively, at the beginning of the time period 
and AV is the change in velocity in time AZ’. Underlined U and Q indicate average values over 
the time AZ’. 

Integration of Eq. 2.10 and substitution of limits yields 


1 1 | . 
The average missile velocity during the time period AZ’ is [V) + (1/2)AV]. Thus the time- 
correction term (see Eq. 2.3) expressed in dimensionless incremental form is 
x 


Eliminating AZ’ between Eqs. 2.11 and 2.12 and solving for AV, the following is obtained 
AV=e+f—Ve? + 2fg + f? (2.13) 


where e = X — Vo 
f= AQ (U ~ Ve)X (A2/U7) 
g=X-U 
The velocity at the end of any step is the summation of the AV’s computed from the be- 


ginning of the integration. 
Incremental distance, AD, was computed by the following: 


AD = [V, + (1/2)AV] AZ’ (2.14) 


where V, refers to the velocity at the beginning of the step. AZ’ is defined in Eq. 2.12. 
Evaluation of acceleration (V = dV/dZ) presented little difficulty since integration was not 

involved. Furthermore, the time-correction term was not necessary because, by definition, 

V is the instantaneous time rate of change in velocity. Thus, the following equation was formed 

from Eq. 2.7: 


(2.15) 


2.3 EVALUATION OF BLAST-WAVE VARIABLES 


2.3.1 General Remarks 


A particular classical blast wave can be completely defined mathematically by the pa- 
rameters of shock strength (p,/py)), duration, and either the velocity of sound or the tempera- 
ture for ambient conditions. Secondary-missile computations were made for selected values 
of shock strength, each of which is applicable to any value of duration (and thus bomb yield) or 
ambient sound velocity between wide limits. 


2.3.2 Dynamic Pressure and Wind Velocity 


Although dynamic pressure (from which wind may be computed) has been measured in ac- 
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tual blast situations, values computed from theoretical considerations were used in this study. 
The reason for this was both the higher reliability and the greater facility for numerical treat- 
ment of the computed parameters over the measured ones. Of the several studies made of the 
blast wave, those made by Harold L. Brode of Rand Corporation»? were found to be most use- 
ful for the present study. The equations’ listed below are empirical relations determined by 
fitting curves to computed blast data. In terminology consistent with the present study, dy- 
namic pressure as a function at shock overpressure and time is given by 


Q=Q, (1 —Z) [Je~%* + Ke~°4] (2.16) 


_ (2.5 P2\ (1 +2107 pt 

where Q= Fre) Geos PE 

J = 1.186 PS for P, < 0.6 

J=1 for 0.6 = P, =1.0 

J = (10* pg4)/(104 + P2) for P, > 1 

K=13 

y = (1/4) + 3.6 P? 

6 =(7+8 PZ) +2 P2/(240 + Pg) 


The relation for wind, or particle, velocity is 
U = Us (1 —2Z) e7 "2 (2.17) 


where U, = (P,)/(1 + PZ) 
y = P4 + 0.0032 PZ 


2.3.3 Overpressure vs. Time and Overpressure Impulse 


Overpressure as a function of time does not enter directly into the computation of 
secondary-missile behavior; nevertheless, it seems appropriate to consider this relation since 
it is the most commonly measured parameter of the real blast wave. Thus, overpressure- 
time can be considered to be the bridge between secondary-missile field data and the com- 
puted data resulting from the present study. 

The following overpressure-time relation was obtained from Brode:'~* 


P =P, (1—T) (ae*! + be") (2.18) 


2.282 (8 + Pg) 


where f= Ce ae a ea A + 0.23 
27.658 + P, + 1.2 PZ + 0.007 PS 
b=l-a 
2 
i= [Pe 15 PE 
1+0.1 Ps 1500 + Pz 
j=9+1.4P, 


Pressure instrumentation used in field work often produces a more accurate measure - 
ment of overpressure impulse than of shock overpressure. Indeed, an improved estimate of 
shock overpressure can be made by making use of the impulse relation described below. 

Overpressure impulse is defined as 


I= es P dt (2.19) 


However, to facilitate integration of Eq. 2.18 in terms of normalized time, the following rela- 
tion was used: 


T =t/tp (2.20) 
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thus, 
dt = t} dT (2,21) 


A combination of Eqs. 2.19 and 2.21 gives 


=tt f' pat 2.22 
I=tpf Pd (2.22) 


Integration of Eq. 2.18 in the manner indicated by Eq. 2.22 yields the following: 


I= Pst Fi (e1 +i -1) +7 (e“J+j- ») (2.23) 


Figure 2.1, a plot of P, in atmospheres as a function of I/t5 also in atmospheres, illus- 
trates this relation graphically. This plot can be thought of as defining the “shape factor” of 
the overpressure-time curve as a function of maximum overpressure. If impulse, I, and dura- 
tion, th, are measured by suitable instrumentation, then a value of shock overpressure can be 
determined from the curve shown in this figure. 


2.3.4 Duration Concepts 


(a) Positive-overpressure Duration. To evaluate the computed motion parameters for a 
particular yield, it is obviously necessary to know the duration of the blast wave (identified by 
peak or shock overpressure) of interest. For this purpose the durations computed from theo- 
retical considerations for free-air conditions, such as those by Brode, are of little value since 
the complex effects of surface reflections are not considered. Thus, the semiempirical rela- 
tions presented in Chap. 3 of The Effects of Nuclear Weapons‘ were used to define overpres- 
sure duration as a function of yield, overpressure, ambient pressure, and the speed of sound. 
Using data for both the surface burst and the “typical air burst,” the following mathematical 
expression was derived 


log th = 5.7995 + (1/3) log W — 0.2957 log pg — 0.0376 log po — log cy (2.24) 


where tp, = duration of positive overpressure in milliseconds 
W = yield in kilotons 
p, = shock overpressure in pounds per square inch 
Pp = ambient pressure in pounds per square inch 
Cy) = velocity of sound in the undisturbed air in feet per second 


The above equation reflects data for the surface burst for shock overpressure (sea-level 
conditions) from 1.68 to 36.7 psi and for the “typical air burst” from 1.86 to 19.7 psi. 


(b) Overpressure vs. Wind Duration. Instrumentation used in past weapons tests was not 
refined enough to establish a relation between overpressure and wind duration. However, the 
theoretical work quoted above has established such a relation. Figure 2.2, derived from 
Brode’s work,'*’ presents the ratio of the wind duration to the pressure duration for overpres- 
pures up to 3 atm (44.1 psi at sea level). It is apparent from this chart that for the higher 
overpressures air-particle inertia has the effect of sustaining positive winds for an appreci- 
able time after the overpressure has become negative. 


2.3.5 Velocity of Propagation of the Pressure Disturbance 


It has been shown that it is necessary to know the velocity of propagation of the pressure 
disturbance in order to compute the motion of objects displaced by blast waves. An easily 
evaluated relation’ that was used for this purpose is: 


x == + V1 + (¥) (2,25) 
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Fig. 2.2-—Ratio of duration of wind to that of overpressure as a function of shock overpressure. 
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Two objections might be raised to the use of the above equation for the purposes of this 
study. First, it applies strictly to the speed of propagation of the shock front, not to pressure 
regions behind the front. Second, it was derived for nondivergent flow; whereas the present 
study applies to divergent flow. In spite of these limitations, the relation was found to be in 
reasonable agreement with work done by Brode!** as quoted in Sec. 2.3.2. 

This, added to the fact that X appears only in the time-correction term (see Eq. 2.7), 
whose effect on the computed value of dV is second order, probably justifies its use in the 
present context. 
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Chapter 3 


COMPUTATIONAL METHOD 


3.1 SCOPE OF COMPUTATIONAL EFFORT 


Because computations were made in terms of dimensionless quantities, it was necessary 
to use only two independent variables: (1) the acceleration-coefficient numeric* (A = apgti/co) 
containing acceleration coefficient, ambient pressure, speed of sound, and duration of positive 
winds (yield dependent) and (2) the shock-overpressure numeric (P, = ps/pp), which also in- 
volves the ambient pressure. Thus, five independent variables, one describing the object dis- 
placed and four describing the blast wave, were effectively reduced to two. It should be pointed 
out that the shock-overpressure numeric (along with the duration variable) represents or de- 
fines three other blast variables that are actually used in the computations; namely, dynamic- 
pressure numeric (Q), wind numeric (U), and propagation-velocity numeric (X), all of which 
are functions of the time numeric (Z = t/t?). 

Thus, for computational purposes, it was necessary to set limits only on A and P,. Con- 
sistent with the scope of the problem stated in Sec. 1.2, the limits arbitrarily set for A were 
0.1 to 9000 and for P, from 0.068 to 1.7 (1 psi to 25 psi for sea-level ambient pressure). Com- 
putations were made for 15 values of P, within the stated range, and associated with each of 
these were 11 values of A (see Table 3.1, columns I and I), making a total of 165 complete 
numerical integrations. 


3.2 GENERAL PLANNING 


Figure 3.1 illustrates some of the considerations in planning for numerical solutions of 
the mathematical model. This plot shows the pertinent blast variables in dimensionless form 
as functions of the time numeric (T = t/th) for a 0.5-atm blast wave. Also shown on this plot 
are velocity (V = v/cy) and displacement (D = d/tic,) computed for an object with an accelera- 
tion coefficient (A = apott/ Cy) of 30. It should be noted that all plotted quantities change most 
rapidly at early times. This means that a stepwise solution should be started using small time 
increments, these to be lengthened as the solution progresses. Also of interest on this chart is 
the U’ curve, which represents the wind numeric at the position of the moving object rather 
than at a fixed position. At T = 0.5, missile velocity was equal to that of the wind, and so the 
integration was stopped. Sixty-two steps were taken to arrive at the solution shown here. 


3.3 STEP SIZE 
As shock overpressures increase, the rate of decay of the blast variables from shock 


values also increases. For mean-value assumptions (i.e., the assumption that the variable is 
constant with a value equal to the mean over a specified time increment) to be equally valid for 


*Numeric is used here to designate a dimensionless quantity. 
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TABLE 3.1—INCREMENTAL VALUES OF THE INDEPENDENT 
VARIABLE AND PARAMETRIC VALUES FOR WHICH COMPUTATIONS WERE PERFORMED 


I I Ii] IV 

P, ty /tt A AT* Tt 
0.068 0.900 0.1 0.0001 0.002 
0.10 0.885 0.3 0.0002 0.004 
0.15 0.875 1 0.0003 0.008 
0.20 0.855 3 0.0004 0.015 
0.25 0.840 10 0.001 0.030 
0.30 0.835 30 0.002 0.060 
0.35 0.805 100 0.003 0.120 
0.40 0.793 300 0.005 0.250 
0.50 0.760 1000 0.007 0.500 
0.60 0.740 3000 0.010 0.750 
0.70 0.720 9000 0.025 1.000 
0.80 0.710 0.050 Final 
1.00 0.675 0.100 
1.30 0.635 
1.70 0.585 


*Ten steps of each AT were used starting with AT = 0,001. See Sec. 3.3 for an explana- 
tion of first four values of AT. 
TT; = times for which computed results were printed out. 


high overpressures, the time increments should be correspondingly decreased. Noting that the 
ratio of overpressure duration to wind duration decreases for increasing overpressures (see 
Table 3.1, column I) suggested the use of a set of time increments in T constant for all solu- 
tions. The AZ values computed for each overpressure (using AZ = tp/t, AT) then decrease as 
desired for the higher overpressure blast waves. 

The first computation in each integration series was made for a time increment, AT, 
equal to 0.001. If the velocity V so computed was greater than 0.1, the solution was discarded; 
then AT values of 0.0001, 0.0002, 0.0003, and 0.0004 were used, in turn, and T = 0.001 was 
arrived at in four steps. Succeeding steps, gradually increasing in size, were then taken until 
the end of the integration (see Table 3.1). If the initial step, AT = 0.001, yielded a velocity less 
than 0.1, the integration proceeded from there without the use of the shorter steps. 

The shortest integration, using the system described above, required 14 steps; this was 
for Ps = 1.7, A = 9000. In general, the number of steps required increased as A decreased; 
e.g., 75 steps were required for A = 3.0 and P, = 0.068 and 81 steps were required for A = 0.1 
and P, = 1.7. 


3.4 MACHINE OUTPUT* 


Since printing out results at the end of each step would have slowed the computation con- 
siderably and also would have produced much more information than could have been utilized, 
it was decided to limit the output of intermediate results to those necessary for the prepara- 
tion of accurate plots. Because of the time-correction term (see Sec. 2.2.2), time at the end 
of any particular step was different for each set of conditions. For simplicity in monitoring 
results and ease of plotting time histories, it was convenient to program output at a set of 
preselected time (see Table 3.1, column IV). Thus, it was necessary to program the com- 
puter to interpolate (linearly) the computed results between time steps to the times selected 
for print-out. 

Special problems arose in the determination of the final values of the computed results; 
i.e., the values occurring at the time when missile velocity and wind velocity were identical. 


*The computer, a CRC-102A, was generously made available by Sandia Corporation. 
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Since missile velocity changes very slowly near the end of the accelerative phase, it was suf- 
ficiently accurate to take final or maximum missile velocity to be that computed for the first 
step where it equaled or exceeded the wind velocity. However, it was necessary to obtain the 
time at which this occurred by interpolating the wind values to a time when they equaled maxi- 
mum missile velocity. Making use of this time, displacement at maximum velocity could then 
be computed. Final acceleration was always, of course, zero. 


3.5 DISCUSSION OF ERROR 


The usable word length of the computer was 36 binary digits or the equivalent of 9 decimal 
digits of input or output. The fixed-point fractional mode of operation required careful scaling 
of all magnitudes (primarily because of the large range of the parameter A) so that sufficient 
significance be retained without the need for rescaling. Binary scaling proved adequately con- 
servative in the attainment of this objective. 

Approximations for square roots and cube roots were obtained with accuracy greater than 
eight decimal places, and the exponential approximation is reported to be accurate to +2 in the 
seventh decimal place. (Cube roots were obtained by the Newton-Raphson method, and the ex- 
ponential, by the rational polynomial approximation.) 

Blast-wave parameters were evaluated from empirical equations that were derived by 
fitting curves to computed data obtained from a blast-wave model.*** Although Brode did not 
make a definite statement regarding the over-all accuracy of the blast model, he did indicate 
some deviation of the empirical equations from the computed data. From this it can be sur- 
mised that computations involved in the present problem were carried out as accurately as 
was warranted by the accuracy of the input blast data as well as by the probity of the missile 
mode! itself (see Sec. 1.2). 

It is noteworthy that computed missile velocity becomes stable by virtue of the number of 
steps involved in each integration; i.e., if for some reason computed missile velocity at the end 
of a particular step is too low, the net wind velocity (U — V) used in the next step will be cor- 
respondingly high, thereby tending to compensate for the original error. As a consequence the 
final solution is not extremely sensitive to the magnitude of the time increments so long as, 
within any particular solution, the steps are sufficiently numerous for the compensatory ef- 
fect to be realized before the computation ends. 
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Chapter 4 


RESULTS: COMPUTED MOTION PARAMETERS FOR 
OBJECTS DISPLACED BY CLASSICAL BLAST WAVES 


The results of the 153 numerical integrations are presented in Table 4.1 in terms of the 
dimensionless parameters. Each integration was made for a specific combination of over- 
pressure (P) and acceleration coefficient (A). Values of missile velocity (V), distance of travel 
(D), and missile acceleration (V) are given for 13 times during the accelerative phase of mis- 
sile displacement. Numbers appearing in parenthesis after V, D, and V are scaling factors 
indicating the number of places the decimal point has been moved to the right. Consider, for 
example, the data tabulated for P = 0.10 and A = 1000 at T = 0.120. In this instance V(6) is read 
as 55677, and thus V = 0.055677. 

At T = 0, the time of arrival of the blast wave, velocity and displacement are zero, but 
acceleration is maximum. “T = Final” is defined as the time after arrival of the blast wave 
when missile velocity is maximum and acceleration is zero. The displacement (D) tabulated 
under “Final” is defined as the total displacement of the object at the instant when the velocity 
is maximum. The actual time when maximum velocity is attained appears in the last column 
under “T final’ 

It should be noted that the time measurements used in this table are normalized with re- 
spect to the duration of the positive pressure phase. Since the duration of positive winds is 
longer than that of the positive pressure, T s;,,; values are sometimes greater than unity. The 
exact value of Ts;,,) is, of course, the time when the missile velocity equals the wind velocity. 

Examples of uses of the tabulated data along with various plots are given in Chap. 5. 
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TABLE 4, 1—COMPUTED MOTION PARAMETERS FOR OBJECTS DISPLACED BY CLASSICAL BLAST WAVES 


Nomenclature (Parameters can be evaluated in any consistent system of units): 


a = sCh im m = mass of moving object tt = positive wind duration 
i 
_ + = i = . 
i = apo ty a P., = maximum overpressure 4 ee os aes 
eee velocity of sound in undisturbed air P, = Pressure of undisturbed air v = velocity of moving object 
eee : n 
Cy = drag coefficient of moving object P= pL [Py V(n} = v r ci x 10 
d = displacement of object S = area presented to wind by object v = acceleration of moving object 
Din) = d/ Cc c,) x 10” tS positive overpressure duration Vin) = ¥ i / co x 10° 
P A T 0 . 002 -004 .008 .015 .030 .060 .120 .250 .500 .750 41,000 Final Trinal 
V (7): 0 88 175 345 635 1218 2253 3912 6333 8820 9941 10296 10312 
. 068 3 D (8): 0 ] 3 l2 44 169 641 2327 8439 25850 47142 700620 78517 1,092 
V (7): 49066 48500 47950 46860 45040 41460 35410 26640 16090 7350 2990 450 0 
V (7): 0 293 582 1149 2109 4038 7433 12801 20435 27852 30821 31449 
10 D7): 0 l 4 14 56 até 766 2750 8303 14963 22561 1,020 
V (6): 16355 16155 15958 15573 14930 13670 11560 8541 4970 2080 693 0 
V (7): 0 877 1741 3433 6278 11930 21674 36500 5616) 72687 77374 77687 
30 D{?7): 0 1 3 12 43 167 625 2219 7761 22595 39619 49734 0,895 
V (6): 49066 48359 47666 46319 44085 39766 32709 22998 12132 3979 645 0 
V (6): 0 291 576 1127 2037 3777 6578 10369 14477 16736 16896 
100. ~D (7): 0 3 10 42 141 537 1952 6613 21514 57390 84057 0,676 
V (5): 16355 15998 15651 14987 13915 11942 8989 5442 2159 293 0 
V (6): 0 864 1692 3245 5674 9907 15714 21879 26410 27340 
300 =D (7): 0 8 31 120 403 1468 4990 15361 44252 94925 0,458 
V (5): 49066 46971 45000 41392 35994 27296 16832 7547 1650 0 
V (6): 0 2776 5253 9478 15143 22942 30579 35809 37497 37500 
1000) D (7): 0 25 98 366 1151 3772 11148 29385 72750 79332 0,270 
V (4): 16355 14546 13015 10583 7675 4336 1807 466 5 0 
V (6): O 7549 13179 20999 28954 36776 41926 43909 43993 
3000 =D (7): 0 71 259 884 2484 7005 17782 41155 56458 0,159 
V (4): 49066 35858 27310 17282 9279 3569 948 75 0 
V (6): 0 17670 26508 35305 4162& 46042 47903 48084 
9000 D (7): 0 175 579 1712 4167 10150 22914 36728 0,092 
V (3): 14720 6565 3687 1618 621 164 z20 0 
V (7): 0 186 370 732 1346 2586 4797 8373 13662 19090 21447 22143 22172 
. 10 3 D (7): 0 1 3 9 35 134 488 1780 5484 10009 14854 16451] 1, 081 
V (6): 30563 10446 10331 10105 9727 8981 7720 5874 3592 1607 617 79 0 
V (7): 0 619 123] 2433 4466 8550 15748 27142 43334 58605 64117 65030 
10) DD (7): 0 l 2 9 30 117 442 1596 5735 17274 30974 44789 0.991 
V (6): 35211 34780 34357 33530 32151 29451 24934 18455 10662 4157 1159 0 
V (6): 0 186 368 125 1324 2507 4527 7550 11421 14366 14980 14990 
30 = D(?): 0 2 7 26 90 345 1288 4545 15714 44931 77637 87959 0.828 
V (5): 10563 10400 10240 9930 9418 8437 6853 4713 2352 630 40 0 
V (6): 0 615 1214 2364 4238 7745 13170 20056 26737 29507 29549 
100 iD (6): 0 ] 2 9 29 109 391 1291 4051 10392 12684 0, 588 
V (5): 35211 34274 33371 31661 28948 24111 17254 9631 3229 166 0 
V (6): 0 1815 3527 6674 11421 19217 28969 38092 43502 44039 
300 D (6): 0 2 6 25 8] 287 940 2759 7546 12672 0, 382 
V (4): 10563 9957 9398 8408 7002 4921 2717 1042 148 0 
Vv (6): 0 5730 10604 18439 28080 39911 49951 55677 56779 
1000 =D (7): 0 52 197 717 2179 £46786 18945 47387 97492 0,220 
V {4): 35211 29769 25479 [£9234 12654 6243 2219 430 0 
V (6): 0 14926 24836 37145 48177 57673 63035 64483 64485 
3000 D (7): 0 140 496 1613 4298 11440 27641 61685 66790 0, 129 
V (3): 10563 6740 4663 2596 1225 408 88 1 0 
V (6): O 32061 44935 56140 63298 67740 69202 69239 
9000 D (7): 0 323 1015 2836 6574 15342 33599 42608 0,075 
V (3): 31690 10438 5105 1957 672 153 9 0 
V (7): 0 137 273 540 994 1916 3574 6298 10428 14752 16649 17259 17290 
215 ] D (7): 0 2 7 26 98 360 1330 4146 7613 11338 13067 41,114 
V (7): 78671 77860 77070 75510 72890 67720 58890 45730 28730 13000 £5130 970 0 
V (7): 0 411 817 1617 2976 5727 10648 18654 30564 42539 47360 48571 48610 
3 D (7): 0 1 6 20 77 294 1072 3927 12109 22033 32571 35280 1,064 
V (6): 23601 23347 23097 22607 21787 20168 17418 13352 8165 3461 1188 93 0 
V (6): 0 137 ate 537 985 1885 3466 5954 9430 12492 13406 13491 
10 = ~D(7): 0 l 5 19, 66 255 963 3469 12410 36978 65553 87180 0,934 
V (6): 78671 77685 76716 74829 71681 65535 55283 40583 22700 7722 1498 0 
V (6): 0 409 812 1595 2903 5459 9741 15930 23339 28069 28614 
30 =D (6): 0 1 6 19 75 276 962 3251 9009 14907 0.737 
V (5): 23601 23188 22784 22005 20728 18311 14505 £9531 4283 787 0 
V (6): 0 1353 =2659 5140 9097 16235 26597 38527 48331 50886 
100 =D (6): 0 1 5 18 62 230 802 2548 7610 18291 0,493 
V (5): 78671 75942 73342 68503 61051 48440 32022 15799 4053 0 
V (6): 0 3959 7601 14071 23308 37233 52609 64838 70129 70248 
300 D (6): 0 4 14 52 167 572 1776 4920 12715 16408 0,310 
V (4): #23601 21681 19979 17107 %413344 8421 £4021 1252 71 0 
V (6): 0 12152 21747 35909 51461 68068 80049 85424 85866 
1000 =—D {6): n 11 41 143 415 1215 3191 7580 11803 0,176 
V (4): 78671 63322 49088 33388 19384 8167 2417 287 0 
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TABLE 4, 1—COMPUTED MOTION PARAMETERS FOR OBJECTS DISPLACED BY CLASSICAL BLAST WAVES (Continued) 


P A 4 0 002 .004 .008 .015 .030 .060 .120 .250 .500 .750 41.000 Final Tr nal 
V (6) 0 29712 46419 64522 78639 89205 94174 94913 

15 3000 D(7) 4) 281 957 2937 7384 18539 42794 78423 0.103 
V (3): 23602 12265 7478 358) 14712 422 67 0 
V (5) 0 5724 7468 8791 9543 9952 10040 

9000 D(7) 4) 594 1768 4656 10315 23180 49210 0.060 
V (3): 70804 15437 6494 2178 671 124 a) 
Vv (8) 0 709 1412 2797 5159 9972 18694 33219 55666 79530 90144 93841 94326 

.20 -3 D(8) 0 1 2 10 34 131 501 1847 6879 21660 39963 59715 75430 1,195 
V (7): 41667 41270 40880 40120 38830 36280 31870 25160 16130 7380 3000 750 a) 
V (7) ) 236 470 932 1718 3318 6211 11005 18345 25994 29258 30281 30336 

1 D (7) 0 1 3 11 44 167 613 2275 7125 13085 19477 22777 1,127 
V (6): 13889 13754 13620 13360 12921 12050 10551 8277 5237 2321 886 169 0 
Vv (7) 0 709 1410 2791 5140 9905 18458 32445 53308 73877 81694 83460 83509 

3 D({7) 0 1 2 10 33 130 497 1818 6679 20600 37382 55103 59218 1,058 
V (6): 41667 41233 40806 39970 38568 35793 31048 23932 14589 5913 1872 108 ) 
V (6) 0 236 469 926 1699 3246 5958 10200 16013 20819 22032 22099 

10 D (6) 4) 1 3 11 43 162 582 2072 6108 10725 13452 0.894 
V (5): 13889 13711 13537 13198 12632 11527 9687 7043 3807 1145 156 0 
V (6) ) 706 1397 2741 4971 9290 16389 26321 37483 43582 43974 

30 =D (6) 4) 1 2 9 33 125 458 1573 5210 14074 20710 0,677 
V (5): 41667 40854 40063 38542 36066 31440 24313 15301 6234 784 0 
V (6) QO 2324 4550 8730 15269 26672 42325 58939 70827 72901 

100 D (6) 0 2 8 31 103 375 1277 3931 11310 22876 0,437 
V (4): 13889 13299 12743 11726 10203 7746 4781 2129 430 a) 
V (6) O 6744 12803 23233 37417 57327 77289 91341 95988 95997 

300 D (6) 4) 6 23 85 268 888 2651 7054 17590 19256 0.270 
V (4): 41667 37336 33631 27665 20384 11787 5046 1336 14 0 
Vv (5) 0 2017 3504 5549 7607 9598 10884 11356 11371 

1000 D (6) ft) 18 66 223 623 1747 4410 10160 13346 0,153 
V (3): 13889 10028 7569 4729 2501 946 245 16 0 
V (5) 0 4669 6954 9200 10795 11887 12336 12373 

3000 D({7) 0 441 1452 4268 10327 25025 56278 87144 0, 089 
V (3): 41667 18053 9993 4311 #1624 423 49 0 
V (5): O 8302 10362 11795 12545 12921 12976 

9000 D (7): ) 867 2491 6329 13659 30061 54150 0.052 
V (2): 12500 1959 745 230 67 10 0 
V (7): 4) 108 215 427 788 1526 2869 5122 8633 12351 13974 14533 14611 

25 .3 D(7): 4) 1 5 20 75 279 1044 3299 6089 93095 11534 1,199 
V (7): 64655 64070 63500 62370 60480 56690 50100 39890 25720 11590 4620 1170 0 
V (7): 0 360 717 1422 2623 5074 9522 16940 28356 40143 45027 46504 46632 

1 D(7?): ) 1 5 17 65 250 925 3445 10810 19837 29489 34777 1.135 
V (6): 21552 21351 21154 20767 20114 18813 16556 13078 8287 3585 1321 244 0 
V (6): 0 108 215 426 784 1513 2823 4971 8168 11245 12346 12566 12571 

3 D({?7) 0 1 4 14 50 195 745 2733 10052 30943 55942 82193 86992 1.045 
V (6): 64655 63999 63352 62086 59957 55734 48467 37432 22617 8734 2533 91 0 
V (6) a) 360 715 1411 2586 4935 9034 15389 23897 30479 31872 31910 

lo =D (6) 0 1 5 17 64 242 866 3061 8918 15513 18373 0.857 
V (5): 21552) 21268 20990 20448 19546 17787 14861 10662 5542 1469 126 0 
V (6) O 1074 2126 4161 7520 13956 24332 38369 53178 60127 60345 

30 D{(6) 0 1 4 14 49 185 674 2284 7414 19566 26056 0. 628 
V (5): 64655 63256 61899 59300 55100 47358 35710 21496 7949 638 0 
V (6) 0 3528 6880 13102 22646 38751 59733 80440 93498 95022 

100 D (6) ) 3 12 45 151 544 1809 5419 15118 26728 0.396 
V (4): 21552 20464 19452 17631 14981 10901 6316 2564 406 ) 
V (5) 0 1015 1906 3393 5326 7874 10230 11727 12110 

300 D(6) i) 9 33 123 383 1231 3559 9182 21572 0,243 
V (4): 64655 56495 49759 39389 27546 14769 5777 1323 0) 
V (5) 0 2959 5005 7645 10121 12330 13642 14043 14045 

1000 ——D (6) tf) 26 94 310 841 2279 5591 12612 14590 0,137 
V (3): 21552 14433 10322 5995 2938 1026 238 6 9) 
V (5) O 6521 9346 11921 13615 14711 15106 15123 

3000 D({7) 9) 617 1974 5614 13207 31215 68967 94176 0. 080 
V (3): 64655 23671 12102 4817 1713 409 32 a) 
V (5) 0 10886 13173 14661 15394 15735 15768 

9000 D({7) O 1144 3200 7928 16810 36491 58110 0.046 
V (2): 19397 2292 807 234 65 8 0 
V (7) 4) 152 302 600 1108 2150 4054 7268 12307 17623 19908 20690 20804 

. 30 .3 © D({7) 0 2 7 27 104 388 1457 4617 8524 12729 16174 1.201 
V (6) 9247 9168 9090 8937 8678 8160 7249 5814 3762 1673 657 170 0 
V (7) ) 506 1008 1998 3690 7146 13443 23996 40294 56964 63694 65668 65841 

1 D({?7) ) 2 7 23 90 347 1283 4797 15073 27633 41028 48756 1, 142 
V (6): 30822 30548 30277 29746 28849 27054 23912 18990 12031 5091 1819 328 0 
V (6) a) 152 302 598 1103 2128 3975 7008 411505 15732 17162 17414 

3 D(6) 0 2 7 27 103 378 1391 4272 7695 11771 1,035 
V (6): 92446 91547 90641 88868 85880 79935 69639 53809 32177 11873 3175 0 
V (6) ) 505 1003 1980 3627 £6911 12616 21374 32820 41136 42616 42633 

10 =D (6) a) 2 7 23 88 332 1186 4162 11986 20680 23309 0.825 
V (5): 30822 30404 29994 29196 27868 25283 20989 14847 7420 1741 84 0 


TABLE 4, 1—COMPUTED MOTION PARAMETERS FOR OBJECTS DISPLACED BY CLASSICAL BLAST WAVES (Continued) 


p A T: 0 .002 .004 .008 .015 ,030 .060 120 .250 .500 8 ,750 1,000 Final Tr nal 
V (6): 0 1507 +2979 £45820 10481 19319 33295 51603 69808 77210 77317 

. 30 30 D (6): 0 l 5 20 67 253 913 3055 9739 25202 30955 0.590 
V (5): 92466 90267 88139 84084 77580 65751 48368 27908 9431 450 0 
V (5): 0 493 959 1812 3097 5199 7811 10227 11592 11700 

100 —D (6): 0 4 16 62 205 726 2369 £6931 18862 30043 0. 366 
V (4): 30822 29024 27372 24446 20304 14203 7778 2916 362 0 
V (5): O 1408 2615 4574 #7018 10069 12707 14243 14554 

300 D {6): rf) 12 45 165 505 1584 4456 11223 23530 0.223 
V (4): 92466 78821 67944 51892 34628 17406 6319 1269 ) 
V (5): 0 4007 6614 £49801 12613 14977 16288 16622 16623 

1000 D{6): 0 35 124 400 1057 2791 £6701 14888 15638 0.125 
V (3): 30822 19210 13089 7141 3298 1080 227 i 0 
V (5): 0 8452 11745 14572 16324 17411 17755 17762 

3000 D (6): 0 80 250 692 1593 3697 8066 10009 0.073 
V (3): 92466 29055 13911 5196 1770 393 19 4) 
V (5): O 13411 15882 17400 18113 18421 18440 

9000 =D (7): O 1523 3982 9527 19826 42498 61565 0. 043 
V (2): 27740 2584 852 236 63 6 0 
V (7): 4) 200 399 792 1465 2845 5380 9689 16489 23665 26730 27787 27954 

235 .3 D(?7): 0 l 3 9 35 135 503 1898 6037 11155 16664 21476 1,214 
¥ (6): 12500 12399 12300 12104 11773 11105 9919 8014 5213 2304 902 244 4) 
V (7): ) 668 1330 2638 4874 9455 17828 31936 53819 76087 84941 87520 87766 

1 D(?): 0 1 2 9 30 117 448 1662 6236 19633 35989 53412 64400 1.156 
V (6): 41667 41314 40965 40280 39119 36786 32661 26086 16555 6909 2433 451 9) 
V (6): 9) 200 399 790 1456 2812 5259 9284 15239 20735 22522 22817 

3 D (6): 9) ] 3 9 35 133 488 1798 5512 9904 15059 1,032 
V (5): 12500 12379 12260 12026 11631 10842 9465 7321 4343 £1549 392 0 
V (6): 0 666 1323 2610 4778 #=+«%9094 16560 27922 42468 52503 54045 54051 

10 D({6): 0 1 2 8 29 114 426 1517 5291 15092 25877 28159 0,802 
V (5): 41667 41088 40520 39416 37577 34002 28070 19607 9475 2010 54 0 
V (6): 0 1986 3922 7647 13727 25146 42895 65500 86877 94593 94626 

30 =D (6): 0 2 6 25 86 322 1156 3827 12011 30607 35439 0,563 
V (4): 12500 12178 11867 11278 10338 8651 6225 3462 1084 29 0 
V (5): 4) 648 1255 2356 3986 6582 9678 12393 13795 13873 

100 -D (6): 0 5 21 79 259 908 2912 8356 22309 32959 0, 345 
V (4): 41667 38931 36446 32114 26131 17658 9221 3240 321 0 
V (5): 0 1836 3376 S812 8748 12247 #%15113 16674 16934 

300 =D (6): 9) 15 57 207 623 1917 £5281 13059 25233 0.210 
V (3): 12500 10412 8801 6512 4174 1991 683 123 0 
V (5): 0 5110 8261 11945 15031 17522 18828 19114 

1000 —D (6): 0 44 153 485 1257 3251 £7681 16550 0,118 
V (3): 41667 24324 15896 8237 3639 1133 220 0 
V (5): 0 10385 14091 17125 18918 19993 20302 20305 

3000 D {6}: 4) 97 298 809 1834 4200 9082 10523 0. 069 
V (2): 12500 3439 1561 555 183 39 l 8) 
V (5): 0 15854 18478 20011 20716 20998 21011 

9000 D (7): 0 1799 4610 10861 22381 47630 64499 0.040 
V (2): 37500 2863 897 242 61 5 0 
V (7): 0 85 170 338 625 1216 2308 4176 7152 10315 11679 12169 12276 

. 40 1 D (8): 0 1 3 ll 38 147 569 2126 8072 25810 47818 71554 99719 1.290 
V (6): 5405 5365 5324 5245 5110 4837 4346 3542 2326 #1035 413 124 0 
V (7): 0 256 5310 1013 1874 3646 6909 12481 21305 30570 34470 35803 36019 

~3 D(?7): 0 } 3 11 44 170 636 2410 7679 14187 21181 27648 1,227 
V (6): 16216 16092 15969 15727 15315 14481 12986 10544 6866 2996 1157 314 0 
V (6): 0 85 170 337 624 1211 2287 44107 46931 9775 10879 #%411191 11220 

1  D (7): 0 l 3 lk 37 147 565 2102 7902 24880 45545 67506 81618 1,159 
V (6): 54054 53613 53177 52320 50864 47924 42680 34193 21638 8843 3026 542 0 
V (6): 0 256 509 1009 1861! 3596 6730 11882 19456 26282 28390 28705 

3 D (6); 0 ] 3 11 44 167 615 2265 6918 12385 18505 1,020 
V (5): 16216 16062 15909 15609 15102 14087 12304 9501 5565 1900 445 0 
V (6): ) 851 1690 3332 6095 11580 21018 35224 52969 64529 66013 66017 

10D (6): 0 I 3 ll 37 143 534 1893 6554 18495 31501 32888 0.776 
V (5): 54054 53278 52516 51035 48573 43790 35884 24677 11475 2164 17 ) 
V (5): 0 254 500 973 1740 3167 £45343 8033 10449 11215 11216 

30 =D (6): 0 2 8 31 107 402 1429 4678 14454 36285 39581 0, 537 
V (4): 16216 15763 15327 14503 13201 10894 7656 4095 1183 13 0 
V (5): 9) 825 1591 2966 4966 8064 11617 14573 15964 16015 

100 —D (6): 0 7 26 99 321 1109 3499 9858 25862 35565 0. 326 
V (4): 54054 50090 46531 40426 32208 21021 10478 3453 261 0 
V (5): O 2318 4220 7151 10567 14457 17496 19044 19255 

300 D (6): 0 19 71 254 753 2271 #4216139 14938 26743 0. 198 
V (3): 16216 13187 10925 7832 4828 2196 713 114 0 
V (5): 0 6310 9999 14135 17440 20014 21289 21529 

1000 D (6): fy) 54 185 576 1465 3724 #8678 17359 O.,111 
V (3): 54054 29529 18538 9154 3887 #£«®41155 207 0 
V (5): 0 12379 16446 19635 21444 22492 22762 22763 

3000 D(6): 0 116 349 930 2080 4709 10108 10979 0, 065 
V (2): 16216 3921 1695 576 185 37 0 


TABLE 4, |—COMPUTED MOTION PARAMETERS FOR OBJECTS DISPLACED BY CLASSICAL BLAST WAVES (Continued) 


Pp A Ts 


. 50 sea 


<<$0< <«0< <U< 
i, a, wr, a, pei, lin, ez. 
ANY AAA KNAW 


10 D (6): 


30 


100 


300 


1000 


3000 


9000 


. 60 


10 


30 


100 


300 


1000 


3000 


~70 ce 


0 


. 002 


. 004 


. 008 


. O15 


. 030 


. 060 


4805 
110 
9858 


14372 
329 
29394 


474} 
109 
9593 


13817 
32] 
27119 


4187 
1002 
7435 


9973 
2558 
13869 


19338 
5716 
1487 


26456 
9130 
801 


30426 
11988 
17 


6264 
140 
13175 


18729 
418 
39247 


6169 
138 
12767 


. 120 


8800 
415 
8189 


26250 
1239 
24275 


8580 
408 
7763 


24382 
E179 
20772 


6863 
3497 
4832 


14261 
8057 
6545 


22955 
15274 
406 


27914 
21309 
84 


. 250 


10842 
1166 
3733 


32238 
3477 
10974 


10422 
1135 
34108 


28792 
3219 
8459 


75199 
9037 
15816 


13977 
18986 
1377 


20176 
32074 
175 


19774 
2014 
7151 


58594 
9992 
20859 


18726 
1941 
6318 


. 500 


15700 
3750 
1652 


46383 
11134 
4744 


14680 
3582 
13611 


38517 
9784 
2720 


89664 
25084 
2500 


28472 
6467 
3065 


83649 
19123 
8669 


26006 
6073 
2366 


.750 


17784 
6962 
660 


52281 
20596 
1823 


16291 
6551 
4542 


41321 
17431 
577 


32092 
11968 
1196 


93707 
35233 
3220 


28589 
11030 
736 


1. G00 


18546 
10427 
209 


54316 
30763 
521 


16739 
9701 
837 


33402 
17885 
388 


97100 
52470 
930 


29257 
16253 
127 


Final 


23493 
67114 
0 


18737 
14825 
0 


54742 
42156 
0 


16786 
11991 
0 


41691 
25665 
0 


91114 
41899 
0 


14736 
46961 
0 


20198 
40074 
0 


23731 
29349 
0 


26158 
18758 
0 


27462 
11778 
0 


28230 
71606 
0 


26325 
20689 
0 


76552 
58498 
0 


23148 
16217 
0 


55989 
32969 
0 


11724 
50288 
0 


18240 
53438 
0 


24247 
43920 
0 


28011 
31570 
0 


30558 
19952 
0 


31917 
12456 
0 


32714 
75464 
0 


33790 
27204 
0 


97878 
74779 
0 


29325 
20284 
0 


Teinal 


0, 038 


1,310 


1,274 


1, 180 


1,010 


0,744 


0, 503 


0. 302 


0. 182 


0, 102 


0. 060 


0. 035 


1, 330 


1,292 


1,182 


0. 988 


0. 709 


0, 472 


0. 282 


0.170 


0. 095 


0. 056 


0. 032 


1, 384 


TABLE 4, 1—COMPUTED MOTION PARAMETERS FOR OBJECTS DISPLACED BY CLASSICAL BLAST WAVES (Continued) 


. 80 


10 


30 


FOO 


300 


1000 


3000 


9000 


10 


30 


100 


300 


1000 


3000 


10 


30 


100 


= 88 8 
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93750 


0 
0 
31250 


aa 


11219 


9593 
106 
41283 


3016 
340 
12228 


7786 
920 
27071 


17452 
2309 
4024 


27050 
4166 
3105 


33526 
6035 
125 


36019 
7126 
30 


17900 
405 
35773 


5352 
1254 
9544 


12411 
3142 


.16841 


23080 
6775 
1650 


30624 
10478 
822 


34651 
13437 
i5 


14462 
656 
13853 


4306 
1959 
40913 


13991 
642 
12921 


39141 
1836 
33452 


10560 
5281 
7106 


20524 
11531 
8281 


30720 
20371 
418 


35976 
27053 
55 


20815 
900 
20826 


6188 
268 
61310 


20000 
876 
19154 


5520. 


2482 
48217 


14384 
6959 
9523 


26600 
14586 
9954 


37976 
24539 
437 


50269 
5392 
14672 


12259 
14375 
2323 


20895 
28030 
1538 


28088 
43749 
_ 4 


24841 
2505 
8970 


7349 
7443 
26071 


23363 
2402 
7802 


61924 
6612 
17594 


14678 
17241 
2594 


24235 
32641 
1524 


35558 
3421 
13312 


10491 
1015 
38463 


33058 
3254 
11289 


8579 
8821 
24295 


19458 
22201 
3219 


30664 
40217 
1610 


65532 69320 


16090 
4164 


14091 
38619 
259 


35480 
7995 
3748 


10399 
23602 
10530 


32088 
7453 
2810 


79544 
19496 
4678 


16607 
45592 
242 


50433 
10863 
5454 


14723 
3197 
15169 


44873 
10002 
3917 


10830 
25618 
6047 


21628 
57530 
244 


.750 


28317 
718 


39802 
14738 
1429 


11590 
43295 
3806 


35065 
13464 
835 


83588 
34073 
714 


56385 
19957 
2063 


16345 
5842 
5416 


48773 
17968 
1133 


11315 
44422 
837 


1,000 Final 


41337 
21965 
458 


11982 
64289 
1075 


35794 
19772 
133 


58504 
29685 
679 


16876 
8654 
1571 


49711 
26300 
185 


69708 
39757 
0 


14200 
57995 
0 


21555 
59489 
0 


28091 
47593 
0 


32087 
33687 
0 


34753 
21070 
0 


36162 
13080 
0 


36985 
79003 
0 


41793 
33526 
0 


12070 
92126 
0 


35862 
24478 
0 


83939 
46544 
0 


16697 
65299 
0 


24822 
65051 
0 


31823 
50915 
0 


36017 
35593 
0 


38782 
22077 
0 


40233 
13640 
0 


41077 
82115 
0 


59207 
47541 
0 


17013 
12941 
0 


49804 
33134 
0 


11352 
59773 
0 


21702 
18883 
0 


31195 
75037 
0 


38982 
56771 
0 


T 


final 


0.978 


0, 690 


0, 454 


0,269 


0. 161 


0, O91 


0. 053 


0, O31 


1, 391 


1,325 


1,185 


0.959 


0, 666 


0,435 


0,256 


0,153 


0, 086 


0. 050 


1, 448 


1,374 


1,203 


0,950 


0, 646 


0.416 


0,243 


TABLE 4, 1—COMPUTED MOTION PARAMETERS FOR OBJECTS DISPLACED BY CLASSICAL BLAST WAVES (Continued) 


P A 


A ie 


0 


- O02 


. 004 


. 008 


2015 


. 030 


. 060 


. 120 


~ 250 


. 200 


, 150 


1, 000 


Final 


7 nai 


V (5): 0 10033 16653 24913 32218 38567 42263 43476 43498 

1, 00 300 6D (6): 0 72 254 827 2198 5839 14111 31574 38927 0.145 
V (3): 93750 59448 40990 22619 10682 3599 846 39 a) 
V (5): 0 22441 31088 38257 42633 45389 46370 46418 

1000 si (6): 0 182 551 1504. 3435 7925 17252 23839 0.081 
V (2): 31250 9401 4355 1595 543 126 10 0 
V (5): O 34854 41589 45220 47022 47849 47935 

3000 D({6): 0 331 857 2038 4227 49043 14623 0. 047 
V (2): 93750 8471 2520 721 194 25 0 
Vv (5) 0 43453 46703 48152 48703 48810 

9000 D(7) O 4754 10898 23754 46671 87711 0. 027 
V(l): 28125 467 109 25 5 4) 
V (7) 0 644 1283 2545 4706 9143 17279 31000 52023 72550 8056! 83400 84435 

l.3 I D(7) ) 2 6 23 89 342 1272 4768 14907 27164 40221 68159 1.523 
V (6): 50904 50500 50100 49312 47968 45231 40280 32103 19851 786! £2931 981 0 
V (6) 0 193 385 763 1410 2735 5153 9199 15302 21098 23258 23963 24163 

-3 1D (6) 4) 2 | 27 102 379 1411 4374 7924 11684 18340 1,435 
V (5): 15271 15144 15017 14769 14345 13486 11940 9413 5699 £2166 761 226 0 
V (6) 0 643 1280 2536 4673 9018 16831 29538 47752 63543 68549 69768 69898 

1 D(6) 0 2 6 23 88 336 1230 4492 13551 24115 35121 45333 1.230 
V (5): 50904 50404 49910 48940 47297 43997 38182 29020 16374 #15421 1540 261 4) 
V (5) 0 193 382 754 1381 2626 4774 #£«°.8022 12117 14960 15538 15578 

3 D (6) a) l 5 19 67 259 970 3446 11967 33980 58322 78039 0.949 
V (4): #15271 15057 14847 14438 13757 12428 10219 7064 3353 781 101 0 
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Chapter 9 
INTERPRETATION OF RESULTS 


5.1 GENERAL REMARKS 


The motion parameters of secondary missiles computed in the present study can be used 
in many different ways. The purpose of this chapter is to point out a few analytical techniques 
that have been found useful. Although the treatment is not exhaustive, such practical subjects 
as weapon scaling, acceleration coefficients, and interpolation techniques are discussed. In 
addition, samples of computed missile data are shown in graphic form. 


5.2 ACCELERATION COEFFICIENTS FOR VARIOUS OBJECTS 


Acceleration coefficient has been defined as the product of the area presented to the wind 
by an object and its drag coefficient divided by its mass. To vivify the meaning of the com- 
puted results, it is necessary to relate values of acceleration coefficient to real OnJects: For 
this purpose Table 5.1 (based on Refs. 1—3) was prepared. 

In regard to the acceleration coefficient for man, it is evident from these data that posi- 
tion with respect to the wind is quite important. Change of position during translation, as well 
as surface-friction effects, serves to complicate any attempt at an exact analysis. With re- 
spect to this problem, it is useful to note the results of an experiment reported by Taborelli 
et al.‘ in which anthropometric dummies, weighing 169 lb dressed and having a height of 5 ft 
9 in., were used in connection with full-scale weapons tests. In one instance a dummy was 
placed on a concrete ramp standing with its back to the oncoming blast wave. The blast pa- 
rameters at this location were: maximum overpressure 5.3 psi, ambient pressure 13.3 psi, 
duration of positive pressure 0.964 sec, and the velocity of sound in the ambient air 1120 
ft/sec. Partial results of the motion picture analysis shown in Fig. 5.1 indicate that the 
dummy was accelerated to 21.4 ft/sec in 0.5 sec after having been displaced about 8 ft. By 
this time the dummy’s position, which was initially vertical, had become horizontal with the 
head toward the oncoming blast wave (see chart at top of Fig. 5.1). Also shown on this figure 
are predicted velocity-time histories for various values of acceleration coefficient. It is in- 
teresting to note that early in the displacement record, the dummy’s velocity corresponded 
closely with that predicted for a man standing broadside to the wind (a = 0.052 ft?/lb, see 
Table 5.1). At later times, because of rotation, the dummy’s increase in velocity with time 
corresponded more closely with that predicted for a prone man aligned with the wind (see 
lower curve on Fig. 5.1). It should be noted that the record obtained for the dummy was ter- 
minated because dust obscured the test area; therefore, the latter part of the record was less 
accurately determined than the initial portion. There was some indication that the dummy’s 
velocity was increasing slightly at the termination of the record. However, if 21.4 ft/sec is 
assumed to have been the maximum velocity attained, then it is possible to determine an ac- 
celeration coefficient that would produce a predicted maximum velocity of the same value 
(21.4 ft/sec) under the same blast conditions. The effective acceleration coefficient so de- 
termined had a value of 0.0268 ft?/Ilb. This:is remarkably close to 0.03 ft?/lb computed - 
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TABLE 5.1—-TYPICAL ACCELERATION COEFFICIENTS (a) 


a, ft/lb Reference 
168-iIb man: 
Standing facing wind 0.052 1 
Standing sidewise to wind 0.022 1 
Crouching facing wind 0.021 1 
Crouching sidewise to wind 0.017 1 
Prone aligned with wind 0.0063 i 
Prone perpendicular to wind 0.022 1 
Average value for tumbling man in 
straight, rigid position 0.030 2 
2i-g mice, maximum presented area 0.38 2 
180-g rats, maximum presented area 0.19 2 
530-g guinea pigs, maximum presented 
area 0.15 2 
2100-¢ rabbits, maximum presented area 0.079 2 
Typical stones: 
0.lg 0.67 2 
1.0 ¢ 0.32 2 
10.0 ¢ | 0.15 2 
Window-giass fragments, '4 in. thick: 
0.1 g, ail orientations 0.78 2 
1.0 g, edgewise and broadside to wind 0.48—0.57 2 
10.0 g, edgewise and broadside to 
wind 0.34—0.72 2 
Steel spheres: 
¥, in. diameter 0.139 8 
¥, in. diameter 0.0696 3 
¥\, in. diameter 0.0398 3 
¥, in. diameter 0.0348 3 
¥,, in. diameter 0.0310 3 


for a tumbling man in a straight, rigid position (see Table 5.1). Using the latter value, a pre- 
dicted maximum velocity of 23.4 ft/sec was obtained in a displacement of 19.5 ft (see plotted 
point in Figure 5.1). The total displacement of the dummy, measured after the event, was 
21.9 ft. This figure included, of course, the distance required for the dummy to come toa 
stop after maximum velocity had been reached, and therefore it cannot be compared directly 
with displacement predicted at the time of maximum velocity. 

Other field studies**® have been made to evaluate the velocities of glass-fragment missiles 
originating from windows facing the oncoming blast wave (the window frames were mounted in 
the open and in houses). It was found in the case of the house-mounted windows that a con- 
siderable portion of the missile sample from each window had velocities higher than could be 
explained if acceleration coefficients noted in Table 5.1 were used in applying the results of 
the present study to the blast situations encountered in the field operations. However, if one 
computed on the basis of a reflected blast wave, the higher missile velocities were satisfac - 
torily explained. Thus, it appears that the somewhat complicated hydrodynamic phenomena 
occurring when a blast wave enters a house by way of a window or windows produces missile 
results equivalent to a shock overpressure more than twice as great as the overpressure ac- 
tually incident upon the house. 

Acceleration coefficients for steel spheres have been included in Table 5.1 for compara- 
tive purposes. For instance, the alphas for a tumbling man and a steel sphere wae in. in di- 
ameter are about the same. Similarly, 4 in. steel spheres and guinea pigs are approximately 
equivalent in so far as translation by blast waves is concerned. Thus, in theory, an “equiva- 
lent” sphere can be found for any irregular object. This concept has been used in weapon- 
effects tests’*® taking advantage of the fact that velocity can be experimentally determined 
more readily for a sphere than for the object it represents. 
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Height above Surface, ft 


Dummy Positions at Various Times 


Velocity, ft/sec 
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Fig. 5.1——Anthropometric dummy translation history, obtained from full-scale weapon test,* com- 
pared with that predicted using various values of acceleration coefficient (see Table 5.1) and the 
computed data in Table 4.1. Numbers adjacent to plotted points indicate measured or computed dis- 


placements. Blast parameters: p_= 5.3 psi = 13.3 psi t” = 0.964 sec, c. = 1120 ft/sec. 
s 3 Py 3 Pp 3 0 
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5.3 WEAPON YIELD AS A BLAST PARAMETER 


Although the motion parameters were evaluated without regard to yield, introduction of the 
latter at this point is both interesting and useful. Employing the weli-known weapon-scaling 
law !?8 applying to given values of shock overpressures, the following relation can be written: 


(to), (ta) W1Po Co 


th _ ty _ ae (Co) (5.1) 
Quantities marked with the subscript “1” are considered to be constant and to have reference 
or “standard” values. The parameters not so marked can take any set of appropriate values. 
The next step is taken from the definitions of dimensionless acceleration coefficient (A) 
and displacement (D). The subscript marking is similar to that for Eq. 5.1. 


ap,tt apott 
Po H _ *Potu (5,2) 
Co /4 Co 
d - ( d } (5.3) 
tico \tuco/; 
Eliminating t{/(t{), between Eq. 5.1 and Eqs. 5.2 and 5.3, in turn, the following is ob- 
tained: 
cosy) | Po |" Gr 
a,=a |—"*| a 2.4 
; | Co (Do) Wi 4) 
W al 
& Po 


Thus, weapon yield replaces duration as parameter. The significance of this transforma- 
tion will be demonstrated in Sec. 5.4. 


0.4 MAXIMUM VELOCITY AND CORRESPONDING DISPLACEMENT 


Data in Table 4.1, along with Eq. 2.24, were used to prepare Fig. 5.2, which shows the 
maximum missile velocity as a function of acceleration coefficient and shock overpressure. 
The tabulated data were made dimensional for a 1-kt burst where the ambient pressure and 
speed of sound were 14.7 psi and 1117 ft/sec, respectively. By use of the transformation equa- 
tion, Eq. 5.4, and the definition of dimensionless velocity, the data on this chart can be made to 
apply to other conditions where W, py, and cy may be different from those used in the construc - 
tion of the chart. 

Consider the translation of a man whose average alpha is 0.03 ft?/lb. For a 1-kt burst 
at a range where the shock overpressure is 1 atm, his maximum velocity is predicted to be 
37 ft/sec (see Fig. 5.2). If the yield were 1000 kt, however, the adjusted alpha, a,, becomes 
0.3. Entering this value for alpha on the same chart, again at the 1-atm curve, a maximum 
velocity of 195 ft/sec is obtained. 

Figure 5.3 shows, for a 1-kt burst, displacement at maximum velocity as a function of 
alpha and shock overpressure, prepared for the same ambient conditions used for Fig. 5.2. 
Continuing the example used above, the man would be displaced 9 ft when maximum velocity 
was reached if the yield were 1 kt. However, if it were 1000 kt, his displacement would be 
28 x 1000 “ = 280 ft (see Eqs. 5.4 and 5.5). 

Figures 5.4 to 5.7 are similar to those described above except that the yields are 20 kt 
and 1000 kt (1 Mt). The charts prepared for 1 kt could be used for all yields, except for limi- 
tations in the range of the abscissa. 
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9.9 ESTIMATION OF MAXIMUM VELOCITY FROM TOTAL DISPLACEMENT 


Estimation of maximum velocity from total displacement can perhaps be described best 
by illustration. Suppose an object of unknown alpha were exposed to blast winds from a 1-Mt 
explosion at a location where the ambient pressure and speed of sound were 14.7 psi and 1117 
ft/sec, respectively, and the shock overpressure was 0.35 atm (5.14 psi). After the explosion 
assume a totai displacement of 100 ft was measured. This measurement includes, of course, 
the distance traveled by the object in stopping after maximum velocity had been reached. By 
use of the plotted data shown in Fig. 5.7, we can only say that the effective alpha must have 
been less than 0.0265. Then, referring to Fig. 5.6, it is determined that for an alpha less than 
0.0265 under these blast conditions, the maximum velocity must have been less than 45 ft/sec. 

The largest source of error in the above estimation, no doubt, embodies the use of total 
displacement as the displacement at the time of maximum velocity, i.e., “overshoot” was neg- 
lected. By means of simple experiments, the amount of overshoot, or stopping distance, could 
be determined as a function of initial velocity. Inclusion of this factor in the estimation pro- 
cedure would result in a smaller but more accurate value for estimated maximum velocity. 

Assuming that the overshoot is known, what are the assumptions involved in estimating 
maximum velocity from displacement? Essentially, it is assumed that the velocity-time his - 
tory of an object is determined by the fact that it traveled a known distance in a known time. 
Mathematically, any number of velocity-time curves could satisfy the known distance-time 
values; however, the most likely one is assumed, in the above procedure, to be of the form of 
computed secondary-missile velocity vs. time for a classical blast wave. It is interesting to 
note that no previous knowledge of alpha is necessary and that an “effective” value is auto- 
matically obtained by the analytical procedure. Effective alpha could be modified in the real 
blast situation by such extraneous influences as ground friction or shielding without seriously 
affecting the accuracy of maximum velocity determined by displacement. However, if the mis- 
sile were in the air at a considerable height at the time of maximum velocity and if the latter 
were fairly high, the estimate of overshoot could be seriously affected. 


5.6 COMPUTED VELOCITY AND DISPLACEMENT FOR PARTICULAR OBJECTS 


9.6.1 Interpolation of Alpha and Overpressure 


Application of the computed motion parameters to specific objects and blast situations 
makes it necessary to interpolate between the values of alpha and/or shock overpressure for 
which computations were made. It has been found that linear interpolation produces results 
sufficiently accurate for most purposes. Graphic interpolation has also been found to produce 
satisfactory results. Charts prepared by such procedures will be presented later in this 
section. 


9.6.2 Velocity and Displacement Predicted for Man and for Glass Fragments 


Figure 5.8 was prepared, primarily, to illustrate the effect of yield on the velocity and 
displacement predicted for a tumbling man and 10-g fragments of glass arising from a window 
in a house (see Table 5.1 @ values). The predictions apply to a shock overpressure of 5.14 psi, 
ambient pressure of 14.7 psi (P = 0.35 atm), and speed of sound of 1117 ft/sec. Also shown on 
Fig. 5.8 are the ranges from Ground Zero as a function of weapon yield where the stated shock 
overpressure could be expected to occur for a surface burst and a “typical” air burst.’ 

For illustrative purposes, the steps taken in the preparation of Fig. 5.8 are outlined. 
Velocity-displacement relations were sought for six yields from 1 kt to 20 Mt. For each yield, 
wind durations were computed for P = 0.35, p, = 14.7 psi, and cy) = 1117 ft/sec, using Eq. 2.24 
and Fig. 2.2. Dimensionless alpha, A, was then computed for each of the yield values used. 
The next step was to obtain velocity-displacement data from Table 4.1 for P = 0.35. For each 
of the times (T) listed in this table, corresponding values of velocity and displacement were 
computed by linear interpolation for the exact value of dimensionless alpha applicable to the 
yield being processed. Then velocity was plotted as a function of displacement for each of the 
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six yields. Velocities were obtained from these plots at given displacement intervals until 
maximum velocity was reached. These data, along with the computed maximum velocities, 
were then used to plot the curves appearing in Fig. 5.8. 

The procedure described above applies strictly to the treatment of the data for man. For 
the 10-g window-glass fragment study, it was assumed that the incident shock overpressure of 
0.35 atm was reflected to 0.80 atm (see Sec. 5.2). 

The plotted data appearing in Fig. 5.8 indicate that the velocity predicted for man is much 
more yield dependent than that for glass fragments. The reason for this is that the window 
glass having a higher alpha reaches wind velocity in a shorter time than does the man and 
therefore utilizes less of the longer duration produced by higher yield. 

It is interesting to note that, for all yields, in only 1 ft of travel the glass fragments have 
already attained velocities greater than 100 ft/sec. Similarly, for all yields greater than 20 kt, 
man is predicted to be propelled at more than 10 ft/sec in just 1 ft of travel. 


9.6.3 Predicted Maximum Velocities and Corresponding Displacements for 1-g Stones 


The purpose of this analysis was to study the interplay of the effects of shock overpres- 
sure and weapon yield on the velocity and displacement of 1-g stones. The results, shown in 
Fig. 5.9, were plotted so that corresponding values of shock overpressure and weapon yield 
could be obtained for the plotted values of maximum velocity and displacement at maximum 
velocity. Data for this chart were scaled from those presented in Figs. 5.2 and 5.3. 

An interesting concept to be derived from this analysis is that of “equivalent” shock over- 
pressures; e.g., a 15-psi blast wave produced by a 1-kt burst is equivalent to a 5.7-psi wave 
from a 20-Mt burst in that both are predicted to propel 1-g stones at maximum velocities of 
200 ft/sec. It should be pointed out, however, that the distances required to achieve maximum 
velocity are quite different. For the 1-kt burst the required distance is about 30 ft; whereas, 
for the 20-Mt burst it is about 300 ft. 

From the data plotted in Fig. 5.9, it can be concluded that both velocity and displacement 
increase with pressure and yield; however, missile velocity is more sensitive to pressure 
‘wind-velocity dependence) than is displacement; and, conversely, displacement is more sensi- 
sive to yield (duration effect) than is velocity. 
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Chapter 6 
DISCUSSION 


Although this study was intended to further the understanding of the secondary and tertiary 
effects of blast on a biological subject, the results are equally applicable to certain other in- 
vestigative efforts, e.g., studies of physical damage resulting from secondary missiles or dis- 
placement. The model that was constructed to describe the motion of objects displaced by the 
blast wave requires no knowledge of the object displaced except its area presented to the wind, 
mass, and drag coefficient, ali of which are assumed to be constant throughout the duration of 
the blast wave. In addition, the model requires that other forces which may be present, such 
as gravity and friction due to the ‘object’s moving over a surface, be negligible in comparison 
with those due to blast winds. 

Thus, like most models, certain simplifying assumptions were made (see Sec. 1.3) in the 
interest of feasibility, simplicity, and uniformity. The determination of whether or not the re- 
sults so computed apply with sufficient accuracy to particular situations is the subject of other 
investigations; in particular, those conducted in connection with full-scale nuclear weapons 
tests, the results of which! will be published soon. These field studies included the measure- 
ment of velocities for stones and spheres placed in open areas and also for fragments of glass 
from windows mounted in houses and in open areas. Experiments were conducted in locations 
where the incident blast wave varied from classical to nonclassical types. One experiment was 
conducted inside a shelter with an open entryway making use of steel spheres to estimate the 
velocities at which man might be translated. Other observations made under full-scale test 
conditions used dogs to assess the hazards of secondary missiles in houses and in open areas’ 
and to evaluate the effects of overpressure and displacement inside protective shelters with 
open entryways. °*4 

The model dealt with in the present study describes the motion of a missile up to the time 
of maximum velocity, i.e., when the missile velocity is the same as that of the wind. For some 
applications a more sophisticated model might be desired which would take into account 
surface-friction forces during both accelerative and decelerative phases of displacement as 
well as the decelerative wind forces that occur after maximum velocity has been reached. 
Such a model would have application to large objects in situations where lofting is not likely to 
occur. Since total displacement, not displacement at maximum velocity, might be predicted, a 
model could be used to interpret field data where only total displacement was measured; i.e., 
total displacement along with appropriate blast data might be sufficient to reconstruct the 
velocity-time history of the object. It is important to note that this technique need not require 
a prior knowledge of the object’s presented area, mass, or drag coefficient. 

Again, it must be pointed out that the missile model described here applies to an ideal or 
classical blast wave. However, it is well known that blast waves may be modified during their 
passage through a building or into a structure®~ and by the properties of the terrain over 
which they pass. 1+!8 Thus, atypical or nonclassical wave forms can and do exist. The impor- 
tant point is that empirical data are at hand for missiles energized by such wave forms. Con- 
struction of a theoretical model to predict the behavior of displaced objects under such cir- 
cumstances can and should be carried out using the experimental data available as a check on 
the analytical procedures. 
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Such a model could supplement the present study of blast displacement of objects and al- 
low extension of such thinking to aid in the estimation of missile and displacement damage to 
man somewhat along the lines of a recent study,’ which tentatively set forth estimated maxi- 
mal ranges for human hazards from missiles and displacement wherein the explosive yields 


were 1 and 10 Mt from an explosive source detonated at the surface of the earth at sea level. 
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Appendix A 


APPROXIMATION METHODS TO SUPPLEMENT 
THE COMPUTED RESULTS 


A.1 GENERAL REMARKS 


The approximation methods previously described required lengthy computations for nu- 
merical solution, although the accuracy attained was satisfactory for a wide range of values of 
acceleration coefficient. It was realized that other approximation methods could be devised 
which would require little computational effort but that these methods would be valid only for 
special conditions, i.e., for short times after the arrival of the blast wave and for very large 
or for very small acceleration coefficients. Results of these approximations would serve a 
twofold purpose, that of extending and that of checking the computed results presented in Table 
4.1. Material presented in this appendix describes the steps taken to accomplish this purpose. 


A.2 EQUATIONS OF MOTION APPLYING FOR SHORT TIMES AFTER ARRIVAL 
OF THE BLAST WAVE | 


For short periods after the arrival of the blast wave, the acceleration an object experi- 
ences may be considered to be constant. Implicit in the above statement is the assumption that 
the blast wind and dynamic pressure do not decay and that missile velocity is small compared 
to both the wind and shock velocities. Thus, stated mathematically in numeric form 


qdv_._ | 
a V=Q,A (A.1) 
which, upon integration from zero to Z time and from zero to V velocity gives 


dD 
V =Q,AZ = az (A.2) 


Integrating again between the same time values and between zero and D distance 
If Eqs. A.2 and A.3 are combined to eliminate Z, the following is obtained 

Vv? = 2Q,AD (A.4) 
Q. was evaluated for P, = 1.7 using the equation presented in Sec. 2.3.2, and Eq. A.4 was used 
to compute V as a function of D for A values of 0.1 and 30. The relations between V and D thus 
computed are shown graphically in Figure A.1 as dashed straight lines labeled A = 0.1 and 


A = 30. The curved solid lines that approach tangentially the above mentioned straight lines 
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were drawn from the computed data presented in Table 4.1. It is interesting to note that at 
T = 0.03, the approximation is still fairly good for the case where A = 0.1 in contrast to that 
for A = 30 at the same time T. 


A.3 EQUATIONS OF MOTION FOR OBJECTS WITH SMALL ACCELERATION 
COEFFICIENTS 


If an object’s acceleration coefficient is sufficiently small, it can be assumed that the 
velocity attained in a blast situation will be small compared to the wind and shock velocities. 


Thus, Eq. 2.7 reduces to 


qv_._ 
a7 = V= AQ (A.5) 


Using zero as the initial value of time and velocity, the above expression can be integrated to 
give 


Zz 
7 =) 
V=A | QdZ == (A.6) 


Equation A.6 can be integrated similarly to give 
D=A {* {[* Qaz az (A.7) 
0 “0 


Combining Eqs. A.6 and A.7 to eliminate A, the following is obtained 


z 
‘7 J, Qadz 
aaa cera rae 
f° f° Qdz az 
0 “0 
The evaluation of the above integrals can be accomplished by use of Eq. 2.16 
: J K 
fe dZ = a.|C e774 -S e762) (Z —1) 
0 
Gee ee eee J K J K 
hich Y cai OL 4452 oP fait 
sa e +—@ ( Ber % ‘| (A.9) 


and 


Zz Z 
J { Q dZ dZ = a,|-( e7 ¥4 +s e~82) (Z —1) 
0 0 
23 _oaz , 2K _6z JK J K 
= ae +s e + (—-+—~--3| Z (A.10) 
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The relations derived above were used to describe V as a function of D for a 1.7-atm blast 
wave. These solutions are shown in Fig. A.1 as dashed straight lines for T values of 0.03, 
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1.00, and 1.71. For comparative purposes the accurately computed data from Table 4.1 were 
used to plot the solid curves, which deviate from the approximation lines at the higher values 
of V and D. It should be noted that for this blast wave (P, = 1.7), T = 1.71 corresponds to 

Z= 1.0, i.e., it is assumed for this approximation that the missile was influenced by the en- 
tire positive phase of the blast wave. This, clearly, is the limiting case since it requires that 
the velocity gained by zero and thus A = 0. In spite of these assumptions, this approximation 
(see dashed line for Z = 1.0, Fig. A.1) is in fair agreement with the “final-value curve” at 
A= 0.1. 

The relation between D,, and A is illustrated in Fig. A.2 for three values of P,. As in 
the previous chart, the solid curves were obtained from the accurately computed data and the 
dashed lines represent approximate relations, which are accurate only for extreme values of 
A. Equation A.7, with limits on Z from zero to 1, was used to compute the approximation lines 
for small values of A. The upper approximation lines appearing on this chart are discussed 
in the next section. 


A.4 APPROXIMATION RELATIONS FOR LARGE ACCELERATION COEFFICIENTS 


It was found by trial that, for missiles with large acceleration coefficients, it was suffi- 
cient to consider the maximum missile velocity equal to the peak wind velocity; however, to 
estimate missile displacement at maximum velocity, it was necessary to compute by approxi- 
mation methods missile velocity as a function of time and to take into account the decay of the 
‘wind with time. 

Wind as a function of time can be approximated (for short times) by a straight-line func- 
tion, thus 


U(t) = Us — BZ (A.11) 


‘where B represents the initial slope of the wind-time curve, being a parameter dependent only 
on Pg. 

Using the basic equation, Eq. 2.7, and the approximation written above, the following is 
obtained 


— BZ =v) Ke 


= Us 
dV =Q,A ( 77 tu, % (A.12) 


In the above equation, Q was approximated by Q,; U by U.,, except when V is subtracted from 
U; and V by U, for the time-expansion term. 

To aid in the solution of the above equation, a new variable is defined, S = —-(U, —- BZ —V), 
whose derivative is dS =dV+BdzZ. After appropriate substitutions, Eq. A.12 becomes 


= =e 5 s?+B (A.13) 


Since it is desired to integrate from V = 0 when Z = 0 to V = U when Z = Z,,, the correspond- 
ing limits were determined for S as follows (see Eq. A.11):; S=—Us when Z=0 and S=0 
when Z = Zy,. Thus, application of these limits to the integrated form of Eq. A.13 yields 


U AQ x 
Zm = fia tan aan (A.14) 
AQ .B*s_ le ial 
X,-Us 
The arc tan factor in the above equation approaches 7/2 as A becomes large. 
The initial slope of the wind-time relation (Eq. 2.17) was found to be 
B=U, (v +1) (A.15) 
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The distance traveled by a missile in reaching maximum velocity is approximately the 
product of Vy (* Ug) and the expanded time (see Sec. 2.2.2) required to reach this velocity. 


Xs 


_ A. 
maxX Us (A.16) 


Dy = UsZ 


Substituting Eqs. A.14 (with the evaluation of the arc tan function indicated above) and 
A.15 into Eq. A.16 yields 


ot Us Xs 
=— = A.17 
Din Us AQs (vy + 1) Ks = Us ( 


where vy is defined as function of P, by Eq. 2.17. 

Equation A.17 was used to plot the approximation lines for large values of A appearing in 
rig. A.2. It is of interest to note that for a given value of A (104, for instance) the approxima- 
tion of D,, is better for strong blast waves (P, = 1.7) than for weak ones (P, = 0.068). This is 
probably because in the strong blast wave the missile gains a higher percentage of the peak 
wind velocity than in the case of the weaker wave. 


A.56 NORMALIZED VELOCITY VS. DISTANCE FOR MISSILES WITH LOW 
ACCELERATION COEFFICIENTS 


A relation that has proved useful is normalized velocity (V/V,,) as a function of normal- 
ized distance (D/D,,). Computed data from Table 4.1 were used to prepare the plots shown in 
Fig. A.3 for three values of A for P, = 0.068 (1 psi at sea level). These plots illustrate that 
the smaller the value of A, the slower is the increase in normalized velocity with increase in 
normalized distance. Hidden in this relation is the fact that, for a given blast wave, missiles 
with the smaller values of A are accelerated over longer times and thus longer distances in 
relation to their velocities. Nevertheless, it would be interesting to determine if there isa 
limiting curve of V/V;, vs. D/Dy for missiles with A values approaching zero. To accom- 
plish this, Eqs. A.6 and A.7 were used in the following manner (remembering that the maxi- 
mum velocity and corresponding displacement are reached at a time, Z, approaching unity if 
the value of A approaches zero): 


(A.18a) 


and 
Zz Z 
p Jo J @42 42 | 
> (A.18b) 
m 
J, J, @aZ az 


The above equations were solved for several corresponding Z values, and the results were 


used to plot the curve of A= 0 in Fig. A.3. 
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Fig. A.8--Normalized velocity vs. normalized displacement for various values of acceleration coeffi- 
cient computed for a shock overpressure of 0.068 atm. Data for the solid curves were obtained from 


Table 4.1, and those for the dashed curve (A = 0) were obtained from approximation methods (see 
Sec. A.5). 
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CIVIL EFFECTS TEST OPERATIONS REPORT SERIES (CEX) 


Through its Division of Biology and Medicine and Civil Effects Test Opera- 
tions Office, the Atomic Energy Commission conducts certain technical tests, 
exercises, surveys, and research directed primarily toward practical applica- 
tions of nuclear effects information and toward encouraging better technical, 
professional, and public understanding and utilization of the vast body of facts 
useful in the design of countermeasures against weapons effects. The activities 
carried out in these studies do not require nuclear detonations. 

A complete listing of all the studies now underway is impossible in the 
space available here. However, the following is a list of all reports available 
from studies that have been completed. All reports listed are available from 
the Office of Technical Services, Department of Commerce, Washington 25, 
D. C., at the prices indicated. 


CEX-57.1 
($0.75) 


CEX-58.1 
($2.75) 


CEX-58.2 
($0.75) 


CEX-58.7 
($0.50) 


CEX-58.8 
($1.00) 


CEX-59.1 
($0.60) 


CEX-59,4 
($1.25) 


CEX-59.13 


($0.50) 


The Radiological Assessment and Recovery of Contaminated 
Areas, Carl F. Miller, September 1960. 


Experimental Evaluation of the Radiation Protection Afforded by 
Residential Structures Against Distributed Sources, J. A. Auxier, 
J. O. Buchanan, C. Eisenhauer, and H. E. Menker, January 1959. 


The Scattering of Thermal Radiation into Open Underground 
Shelters, T. P. Davis, N. D. Miller, T. S. Ely, J. A. Basso, and 
H. E. Pearse, October 1959. 


AEC Group Shelter, AEC Facilities Division, Holmes & Narver, 
Inc., June 1960. 


Comparative Nuclear Effects of Biomedical Interest, Clayton S. 
White, I. Gerald Bowen, Donald R. Richmond, and Robert L. 
Corsbie, January 1961. 


An Experimental Evaluation of the Radiation Protection Afforded 
by a Large Modern Concrete Office Building, J. F. Batter, Jr., 
A. L. Kaplan, and E. T. Clarke, January 1960. 


Aerial Radiological Monitoring System, I, Theoretical Analysis, 
Design, and Operation of a Revised System, R. F. Merian, 
J. G, Lackey, and J, E, Hand, February 1961. 


Experimental Evaluation of the Radiation Protection Afforded by 
Typical Oak Ridge Homes Against Distributed Sources, T. D. 
Strickler and J. A. Auxier, April 1960. 


